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We calculate the critical temperature of a superfluid phase transition in a polarized Fermi gas 
of dipolar particles. In this case the order parameter is anisotropic and has a nontrivial energy 
dependence. Cooper pairs do not have a definite value of the angular momentum and are coher- 
ent superpositions of all odd angular momenta. Our results describe prospects for achieving the 
superfluid transition in single-component gases of fermionic polar molecules. 



I. INTRODUCTION 



The recent success in observing quantum degeneracy in ultra-cold atomic Fermi gases |^-^] stimulates a search for 
gaseous Fermi systems with an achievable temperature of superfluid phase transition. The ideas based on Cooper 
pairing for a short-range Van der Waals interaction between atoms ||5[ require a simultaneous trapping of at least 
two different fermionic species, with a rather severe constraint on their relative concentrations. The situation is 
different for Fermi gases of dipolar particles. Being electrically polarized, these particles interact via long-range 
anisotropic (partially attractive) dipole-dipole forces, fn the ultra-cold limit, the dipole-dipole scattering amplitude 
is energy independent for any orbital angular momenta This opens prospects to achieve the superfluid pairing 
in a single- component Fermi gas, where only scattering with odd orbital momenta (negligible in the case of Van der 
Waals interactions) is present. These prospects are especially interesting as in single-component fermionic gases the 
Pauli exclusion principle provides a strong suppression of inelastic coUisional rates (see Q). Hence, one can think of 
achieving significantly higher densities than in Bose gases. 

Possible realizations of dipolar Fermi gases include an electrically polarized gas of polar molecules as they have 
large permanent electric dipoles. The creation of cold clouds of polar molecules has been recently demonstrated in 
experiments with buffer-gas cooling ||ll|] and in experiments based on deceleration and cooling of polar molecules by 
time-dependent electric fields . Another option is to create a gas of atoms with electric dipole moments induced 
by a high dc electric field or by laser coupling of the atomic ground state to an electrically polarized Rydberg state 

. 

The p-wave Cooper pairing in a polarized dipolar Fermi gas has been discussed in ||14| , |10|| , and the corresponding 

critical temperature has been estimated by using the standard BCS approach. In this paper we calculate the exact 
value of the critical temperature and find the energy and angular dependence of the order parameter. For this 
purpose we consider the Cooper pairing for all possible scattering channels. These channels are coupled to each 
other by the dipole-dipole interaction, and the Cooper pairs prove to be coherent superpositions of contributions of 
all odd angular momenta. In order to find the prc-cxponential factor for the critical temperature, we perform the 
calculations to second order in perturbation theory along the lines of the approach of Gor'kov and Melik-Barkhudarov 
(CM approach) 



II. GENERAL EQUATIONS 



We consider a spatially homogeneous single-component gas of fermions having a dipole moment d oriented along 
the z-axis. The Hamiltonian of the system has the form 



*LKB is an unite de recherche de I'Ecole Normale Superieure et de I'Universite Pierre et Marie Curie, associee an CNRS. 
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where iI]{y) is the field operator for fermions, Vd(r) = {d^/r^){l — 3cos^0r) the dipole-dipole interaction, Or is the 
angle between the interparticle distance r and the ^-axis, and /i is the chemical potential. In Eq.(^ we omit the 
contribution of the p-wave scattering due to the short-range part of interparticle interaction, since this contribution 
is small in the dilute ultra-cold limit. 

For a single-component Fermi gas, the Cooper pairing is possible only in the states with an odd angular momentum 
I. On the other hand, the anisotropic character of the dipole-dipole interaction leads to coupling between Cooper 
pairs with different values of the angular momentum. Therefore, the problem of supcrfluid pairing requires us to 
consider states with any odd I. 

The critical temperature Tc of the supcrfluid transition and the order parameter A can be found from the gap 
equation in the momentum representation [ p^ , pT[ : 



A(P) 



dp' 

{2TTh)'' 



tanh(£;(p')/2r) 



A(p') 



(2) 



Here -E(p) = A2(p) -I- [p^ /2m — /i)^, and we assume the order parameter to be real. The function V^(p, p') = 
Vd(p — p') + (5y(p,p'), where Vd(q) is the Fourier transform of the dipole-dipole interparticle interaction potential 
Vd{v): 



An 



d^{icos^{e^) - 1), 



(3) 



with 0q being the angle between the momentum q and the z-axis. The quantity 5T^(p, p') originates from many-body 
effects and is a correction to the bare interparticle interaction Vd- The leading corrections are second order in Vd and 
the corresponding diagrams are shown in Fig.l (see Ref. They describe the processes in which one of the two 

colliding particles polarizes the medium by virtually creating a particle-hole pair. In la the particle-hole pair then 
annihilates due to the interaction with the other colliding particle. In lb, Ic and Id the hole annihilates together 
with one of the colliding particles. In lb and Ic the particle-hole pair is created due to the interaction of the medium 
with one of the colliding particles, and the hole annihilates with the other colliding partner. In Id these creation and 
annihilation processes involve one and the same colliding particle. 
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FIG. 1. The lowest order many-body corrections to the effective interparticle interaction. 

For temperatures just below Tc the order parameter is small and the gap equation is equivalent to the Ginzburg- 
Landau equation for the spatially homogeneous order parameter. This equation can be obtained by expanding the 
rhs of Eq.(0) in powers of the order parameter A(p): 
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where K{p) = tanh(f /2r)/2^, and ^ = p'^/2m - /i. 

The occurrence of the Cooper pairing is associated with the existence of a nontrivial solution of Eq. (Q) for temper- 
atures T < Tc- In order to find the value of the critical temperature one can neglect the second, nonlinear term in the 
square brackets in the rhs of Eq.(^) because for T — > Tc the order parameter A ^ 0. The corresponding linearized 
gap equation also provides us with the momentum dependence of the order parameter, whereas the nonlinear term 
determines the absolute (temperature dependent) value of A. 

The integral in Eq.(^ diverges at large momenta. The divergency can be eliminated by expressing the bare 
interaction Vd in terms of the vertex function ( scattering off-shell amplitude) r(i?,p,p'). This is similar to the 
well-known procedure of renormalization of the scattering length in dilute gases of Bose or Fermi particles interacting 
via short-range forces One may choose any value of E, and for simplifying our calculations we select E = 0. 

Then the vertex function r(0,p,p') = Tii{p,p') obeys the equation 

Taip, p') = Vd{p -P')- J J^^^dip, q)Ko{q)Vd{q - p'), (5) 

with Ko{q) ~ raj . We will confine ourselves to the second order in perturbation theory with respect to Vd- Omitting 
higher order corrections, the rcnormalizcd linearized gap equation reads 

= " / (^r,(p,p'){W)-^o(p')}A(p')- 1 ^^<5nP,P'W)A(p'). (6) 

In the dilute ultra-cold limit only small momenta p and p' are important. We thus have to find the scattering 
amplitude for ultra-cold particles, in the presence of the dipole-dipole interaction between them. 



III. SCATTERING AMPLITUDE IN THE ULTRA-COLD LIMIT 



The anisotropic and long-range character of the dipole-dipole interaction {Vd oc 1/?''^) ensures that in the ultra- 
cold limit all partial waves give an energy independent contribution to the scattering amplitude For any orbital 
angular momentum / one has ~ ~ ATih^r^/m, where the quantity r, ~ md? /Ti^ plays a role of the characteristic 
radius of interaction for the dipole-dipole potential. For the interparticle separation r ^ the potential Vd(r) does 
not influence the wave function of the relative motion of two colliding particles and the motion becomes free. The 
ultra-cold limit requires particle momenta satisfying the condition 

pr^/h < 1. (7) 

The anisotropy of Vd directly couples scattering channels with angular momenta I and /±2. Thus, strictly speaking, 
all even-Z (odd-/) channels are coupled to each other, whereas the scattering with odd angular momenta remains 
decoupled from that with even momenta. 

There are two contributions to the scattering amplitude. The long-range contribution comes from distances r > fi/p 
and gives Td ~ cP for all angular momenta in the incoming and outgoing channels, allowed by the selection rules. 
This contribution can be found by using the Born approximation. The short-range contribution comes from distances 
r < r*. For the scattering with even I, due to the presence of the s-wave, we have again F^ ~ (P or even somewhat 
larger because of the so-called shape resonances ||2^. Under the condition (|^), the contribution of the s-wave to the 
wave function of the relative motion at distances r is independent of p. This leads to an energy independent F. 
However, it depends on a detailed behavior of the interparticle potential at short interparticle distances. Thus, for 
even / one can not make a general statement on the value of F. 

In the case of identical fermions only odd orbital angular momenta are present. Then the short-range contribution 
is much smaller than the long-range one. We will demonstrate this for the p-wave on-shell scattering amplitude, 
omitting the coupling to the channels with other odd I. For / = 1 and nii = in both incoming and outgoing 
scattering channels (m; is the projection of I on the z axis), the dipole-dipole potential Vd(r) averaged over the angle 
Or is equal to Vd = — 4d^/5r^. In order to analyze the short- and long-range contributions to the scattering amplitude, 
we consider the relative motion of particles in a truncated potential V{r) = Vd{r) for r < rp, and V{r) — for r > tq. 

The truncation radius tq is selected such that tq ^ The Schrodinger equation for the wave function of 
the relative motion reads 

f (P 2d 6 \ ,,,,,,,,,, p2 , 
m \ dr^ r dr J m 
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At distances r h/p we may put p = in Eq.(^). Then for r < tq we use the weh-known procedure of finding an 
analytical solution for ip{r) in power law potentials [pl| . Assuming r 3> r*, this gives 



tp{r) oc ^— + const^ , 



(9) 



where the constant term is independent of r^. At r > rg the motion is free and ipir) depends explicitly on the 
scattering phase d. The solution, which for r — > oo takes the required asymptotic form (h/pr) cos {pr /h + S)^ at 
r <ti h/p becomes 

Equalizing the logarithmic derivatives of the wave functions and (^0|), we immediately obtain S ~ p'^r^r^/Ti^ and 
find that tp{r) oc p for r ^ h/p. The scattering amplitude then proves to be = / il){r)V{r)d^r ^ (P{pro/K)'^. The 
short-range contribution to F, that is the contribution from distances r ^ r*, is obtained from this relation by simply 
putting ro '--^ r,. 

We then increase tq and make it much larger than fi/p. At distances r ^ h/p ^ r^. the potential V is much smaller 
than the kinetic energy term in the Ihs of Eq. (^) . For the contribution of these distances to the scattering amplitude 
the Born approximation gives Tfi ^ cP. We thus see that the short-range contribution to the scattering amplitude is 
small as (jpri,/hY compared to the long-range contribution coming from distances r ~ h/p. 

This has two important consequences. First, a detailed shape of the interaction potential is not important for the 
scattering amplitude as the latter is determined by the long-range contribution. This contribution is obtained in the 
Born approximation and depends only on the value of the dipole moment. Second, we may include the second order 
Born correction to the amplitude. This correction is of the order of cP{pr^:/h) and still greatly exceeds the short-range 
contribution. 

In the second order Born approximation for the off-shell scattering amplitude Fc;(p, p') we have 

Frf(p,p')-Ki(p-p')-y ■j^^Vdip-q)Ko{q)Vdiq-p'), (11) 

where the first and second terms in the rhs of Eq.(pT|) are first and second order in Vd, respectively. The integral for 
the second order correction to the scattering amplitude in Eq.(pl1) is formally divergent at large q. This is the same 
non-physical divergency as in the case of short-range interactions ||l8| , p^ , and it will be eliminated in the calculations 
of the order parameter and critical temperature (see Section V) . 



IV. CRITICAL TEMPERATURE IN THE BCS APPROACH 



In a quantum degenerate Fermi gas characteristic momenta of colliding particles are of the order of the Fermi 
momentum pp = {6Tr'^n)^^'^h {n is the gas density). Then, with r* ^ iruP /h^, the condition (^ of the ultra-cold limit 
for interparticle collisions can be written as 

nd'^/ep < 1, (12) 

where ep = P% /2m is the Fermi energy. The Ihs of Eq. (|l2|) is the ratio of the mean dipole-dipole interaction energy 
(per particle) to the Fermi energy. As in the case of short-range interactions [^,|9| , this is a small parameter of the 
many-body theory. It is the condition ( |l^ that allows us to omit the contribution of higher order diagrams and use 
the renormalized gap equation (^. 

Generally, in dilute Fermi gases the critical temperature is exponentially small compared to the Fermi energy ep- 
The exponent is inversely proportional to the Fermi momentum pp and is determined by first order terms in Vd- The 
account of the second order terms provides us with the pre-exponential factor. 

We first calculate A(p) to first order in Vd and find the correct exponent in the dependence of the critical temperature 
on the particle density. For this purpose we should keep in Eq.(H) only the terms which are first order in Vd- This is 
the first term in the rhs of this equation, with Td{p, p') = Vd{p — p')- Then, Eq.(||) can be rewritten in the form 

f°° f dn' 

A(^, n) = - / de'(tanh(r/2r)/2n / —R{^, n; e, n')A(e', n'). (13) 
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Here n = p/p, and 



R{t n; n') = z.(e')rrf(p(C)n,p(C')n') (1 - C'/iC' + tanh(e72T)) , 

where — mp(^)/27r^?i^ is the density of states at energy ^ + fi. The chemical potential fi is equal to the Fermi 
energy: n^ep- 

The main contribution to the pairing comes from the states near the Fermi surface, where |^| , |^'| <C sf- In order to 
single out this contribution in the rhs of Eq.(|l3|), we introduce a characteristic energy cU that obeys the constraint T <C 
cJ, and is of the order of the Fermi energy. We then divide the integral over ^' in Eq.(|lJ) into two parts: (a) the integra- 
tion of n; 0, n') A(0, n') from -tZJ to ZU, and (b) the integration of n; n')A(^', n') - n; 0, n')A(0, n')] 
from —u) to ZU, plus the integration of n; n')A(^', n') from —ep to — tU and from tU to oo. In part (a) we use 
the asymptotic formula 



dC'(tanh(^72r)/2^') w In 



2 exp(7)aj 



where 7 = 0.5772 is the Euler constant. In part (b) we replace tanh(^'/2T) by the step function (omitting the 
unimportant contribution from a narrow interval |^'| ^ T ^ tU) and integrate in parts. As a result, Eq. (p^) takes the 
form 



A(^,n)--ln 



2 exp(7)a; 



i?(e,n;0,n')A(0,n') + i d£,'\n^^^^j ^ n; n')A(^', n')} , 

(14) 



where the first and second terms in the rhs come from parts (a) and (b), respectively. 

One can easily see that the second term in Eq.(|lJ) is small as 1/ ln(2 exp(7)uJ/7rT) compared to the first one. 
Therefore, the second term is only important for the pre-exponential factor in the expression for the critical temper- 
ature and will be omitted in this Section. This is equivalent to the commonly used BCS approach where the kernel 
i?(^, n; n') is replaced by i?(0, n; 0, n') for |^| , |^'| < ZZJ and by zero otherwise. Putting ^ = in Eq.(p^, we obtain 
the following equation for finding the critical temperature: 



A(0, n) = - In 



2 exp(7)a; 



j ^i?(0,n;0,n')A(0,n') 



(15) 



The anisotropic character of the scattering amplitude leads to a nontrivial angular dependence of the order parameter 
A(0, n). In order to analyze the possibility of pairing we expand A(0, n) in terms of a complete set of eigenfunctions 
(/is(n) of the integral operator with the kernel i?(0, n; 0, n'): 

oc 

A(O,n) = ^A,0,(n), (16) 

^i?(0,n;0,n>,(n') = A,(/.,(n), s = 0,l,... (17) 

The functions 0s (n) are normalized by the condition J ((in/47r)(/)^(n) = 1, and they are labeled by the index s in such 
a way that the eigenvalues A.; < A^+i. Then Eq.(p^) reduces to a set of equations 

A. (l + A.ln^^^)=0. 

The appearance of a nontrivial solution for A(0, n) below a certain critical temperature requires the presence of at 
least one negative eigenvalue A^. For a single eigenvalue As» < 0, the critical temperature immediately follows from 
the condition (1 -I- As* ln(2 exp(7)ZI7/7rTc)) = 0, and we have A^* 7^ and Ag = for s 7^ s* . In the case of several 
negative eigenvalues Aj < 0, one has to choose the solution that corresponds to the lowest eigenvalue as it gives the 
highest critical temperature. 

Using Eq.(^) one finds that negative As correspond to eigenfunctions (/)s which are independent of the azimuthal 
angle (p. This means that only spherical harmonics with zero projection m of the angular momentum I appear in their 
decomposition. For these functions the kernel i?(0,n;0, n') can be reduced to its average over the azimuthal angles 
ip and Lp' . Using Eq.(||) for Td{pF^TPF^'), we obtain 
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i?(O,cos6l;O,cos0') = 27r— f - |cos g - cos 6i'| - ij , (18) 

where 9 and 9' are the polar angles for the vectors n and n', and n is the gas density. Note that the first multiple 
in the rhs of Eq.(|l8|) is a small parameter of the theory, given by Eq.(|l^ and representing the ratio of the mean 
dipole-dipole interaction energy to the Fermi energy ep- 

Keeping in mind that due to the Pauli principle only odd angular momenta are present, we obtain the solutions of 
Eq.(g|): 

0,(n) = V2sin(-(l + 2s)cos(f?)j , \,^- — -^——^. (19) 

The lowest eigenvalue is Aq = —12n(P /ttef- Therefore, the angular dependence of the order parameter will be 
characterized by the function 4>o{n) (see Section VI for details). The critical temperature is then given by 

2exp(7)^ / 1 \ 

In the BCS approach the pre-exponential factor (cU) remains undetermined. One can only argue that it is of the 
order £p. We thus have 



In Ref. [|lO| the exponent in the expression for T^'~"^ is only expressed in terms of the scattering amplitude which 
should be found from the solution of a set of coupled equations. The estimate for this exponent in Ref. |1^ takes into 
account only the p — p scattering channel and contains a numerical error. 

In order to find the pre-exponential factor one has to include the contribution from the second term in Eq.(^), 
together with the second order corrections to the eigenvalue Aq. These corrections originate from the second order 
many-body effects and from the second order corrections to the scattering amplitude, described by the second terms 
in Eqs.(||) and ( pi] ) respectively. 



V. GM APPROACH. THE CALCULATION OF THE PRE-EXPONENTIAL FACTOR 



We now proceed with the calculation of the pre-exponential factor in the expression (|20|) for the critical temperature. 
We first consider the contribution of the second term in the rhs of Eq. (p^) , which is logarithmically small compared 
to the already calculated first term. For this purpose we specify the value of U by the condition 



Then, using Eqs.(^ and (0) we obtain the following expression for lu: 



This definition of ZO immediately leads to Eq. (|2^) for the critical temperature and allows us to rewrite Eq. (Q) in the 
form 

A(^,n) = yj ^i?(C,n;0,n')A(0,n') 

f dn' 1 f°° d^' ( „ i?(e, n; 0, n') A(0, n') f dm /" dna , , , , ^ ,1 

(24) 



(22) 
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where the second term in the rhs is proportional to the small parameter of the theory nd? /ep and can thus be 
considered as a perturbation. This follows from the fact that the bracket in this term vanishes for ^' ^0. As a result, 
in contrast to the first term of the rhs, the second term does not contain the large logarithm ln(a7/r) ~ Aq^ ~ {ep/rufi). 

The leading contribution to the angular dependence of the order parameter on the Fermi surface comes from the 
term with s = in Eq.(|l6|): A(0, n) — Ao(/)o(n). Therefore, to the leading order in ndP /ep, the solution of Eq.(24) is 

A(e, n) « ^ / ^^(C, n; 0, n')Ao0o(n'). (25) 
After substituting this expression into Eq.(^) and performing a numerical integration, we obtain 



= 0A2ep, (26) 



Corrections to Eqs.(p5|) and ( p6D are related to the terms with s 7^ in Eq.(|l6|), and from Eq.(p4[) we find that the 
quantities As^o ~ ^o{n(P /ep). These corrections lead to the relative contribution of the order of n(P /ep to the 
pre-exponcntial factor for the critical temperature and hence will be neglected. 

We now calculate the contributions from the second terms in Eqs.(^) and (pi]). As one can see from Eq.(p^), these 
terms result in the correction for the eigenvalue Aq: 



'^Ao^KO) '^U^)^6V{ppn,ppr^)- j ^-^F^p - q)7^o(g)T^d(q - P')} 0o(n'). 



(27) 



The first term in the integrand of Eq.(p7|) originates from many-body effects, and the quantity (5F(p, p') is shown 
diagrammatically in Fig.l. The analytical expressions for the diagrams in Figures from la to Id, read: 

,^„(p,p.,./jji, "0 + p-/^)-f'('--p-/^) ,,^(p_). 



(27r)3 ?q+p_/2 - ^q-p_/ 



2 



KAK '\ I iV(q + P-/2)-iV(q-p_/2) ^^^ 
5V6(p,p) = - / -— T Vd(p_)yrf(q-p+/2), 



(27r)3 




^q+p_/2 


- Cq-p_/2 


dq 


iV(q 


+ P-/2) 


-iV(q-p_/2) 


(27r)3 




Cq+p_/2 


- ^q-p-/2 


dq 


iV(q 


+ P+/2) 


-7V(q-p+/2) 


(27r)3 




'^q+p+/2 


^ ?q-p+/2 



^^c(p,p')-- / y^'^-' ' V ,(p_)Frf(q + p+/2), 

J (^■^y Cq+p_/2 -4q-p_/2 

f-.rf f dq iV(q + p+/2)-jV(q-p+/2) ^^^ 

5Vd{p,p) = - / 7— ^yd(q- P-/2)Vrf(q + p_/2). 



Here p± — p±p', and A^(p) is the Fermi-Dirac distribution at zero temperature. The integrals related to the first term 
in the rhs of Eq. (p7|) , with 6V{p,p') — J2a=a d '^^"(P' P')i were calculated numerically by using the Monte-Carlo 

method. Each of the terms SVa{p,p') provides a correction SX'^^ — [iy{0)And'^ /3]'^_na- For the coefficients rja we find 
7]a = 0.19, rjb = ?/c = —0.08, and rjd = 0.42. Thus, the first term in the rhs of Eq.(|27|) gives the correction 



<5A«=0.45(K0)^"^^ 



,2^ 2 



The second term in the integrand of Eq. ( p7| ) comes from the second order correction to the scattering amplitude 
Td- For the correction to Aq, originating from this term, our numerical calculation gives 

^A^^^=-0.86(k0)^)'. 

Note that the function (j)o{n) is odd with respect to cos 9. For this reason, the integration over dn and dn' eliminates 
the formal divergency of the integral over dq at large q. 
The total correction to the eigenvalue Aq is then 

^0 = + ,5A('^ = -0.41 (iy(Q)^) . (28) 



3 

On the basis of Eqs.([20|), ( p6| ) and (|2^), we obtain the final expression for the critical temperature: 

Tc = ^^H^Eill X 0.42ei.exp(-l/|Ao + (5Ao|) ~ lA4.epeyiY>{~iiep/\2nd^). (29) 

TT 

It is worth noting that if we include only the p — p scattering channel the exponent in Eq. (|2^) will be larger by a 
factor of IO/tt^. The pre-exponential factor becomes then larger by a factor of 1.1. 
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VI. ANISOTROPIC ORDER PARAMETER 



In order to find the temperature dependence of the order parameter for T < Tc, we have to include the nonhnear 
term in the gap equation (|l5|). This term can be written as 



dp' 
(2^ 



■Vd{p - P') 



tanh(C72T)' 



cosh2(e/2r) e/2r 



where we neglect the many-body correction to the interparticle interaction. The expression in the square brackets 
vanishes as |^'| for |^'| —> oo. Therefore, the main contribution to the integral comes from the region of small 
i.e. from p' close to the Fermi momentum pp. This allows us to write 



1 



cosh"(e72r) e/2T 



2^ 



/ ^V,ip-pFn')A'{ppn') J 



tanh(^V2T) 



a^(pO 

tanh(^72r) 



cosh2(e/2r) 



87r2T2 

where C(z) is the Riemann zeta function. As a result, to first order in Vd the nonlinear gap equation reads: 

27CJ /-dn'^^^ 7C(3) ^ ^n' ^_ ^ 



A(e,n)^-ln-L,/_i.(^,„;0,n')A(0,n') ^^^^ 



— i?(e,n;0,n')A3(0,n') 



^A, (l + A,ln|^) ^,(n) 



s=0 



7C(3) 



Y^XsMn) I ^C,%,,3A,,A,,A,3 I =0. 

{S^} 



s=0 



The coefficients Cf ^ ^ follow from the relation 



S 

For temperatures below Tc, satisfying the inequality (Tc — T)/Tc <^ 1, Eq.(^l|) can be rewritten as 



(30) 



With the order parameter A(0, n) from Eq.(|T|), where now A^ = As(T) and A^(r) for T ^ T^, Eq.@ takes 
the form 



(31) 



s=0 



Ao — As Tc — T 
+ A 



Ao 



Tc 



0.(n) 



7C(3) 



^ s— 



^A,0,(n) I ^C:^,^,3A,,A,,A,3 | = 0, 



(32) 



where we neglect higher powers of {Tc — T)/Tc- It can be easily seen from Eq.(^2|) that for T ^ Tc one has 
Aq '-^ {Tc — T)^/^, and that A^ with s > are either equal to zero or proportional to {Tc — T)^/^. Therefore, the 
equation for the leading coefficient Aq is: 

Tr-T 7C(3) 3_ 



Tc 



ilT^T^ 



where the coefficient Cqqq is equal to 3/2. We thus obtain the following expression for the order parameter on the 
Fermi surface (f = 0): 



v/2TC(3) 



A(0,n) = _i^T,J^V^0o(n) = 2.5 r,J^V^0o(n); « 1. 



Tc 



Tc 



Tc 



(33) 



For ^ 7^ 0, i.e. p pp, the order parameter can be calculated by using Eq.(p5|). Fig. 2 shows the numerically 
calculated dependence of the order parameter on the modulus of the momentum p for various values of the angle 
9 between the vector p and the direction of dipoles. Note that for both s and p-wave pairing due to a short-range 
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interaction, the order parameter is momentum independent for p satisfying the condition of the ultra-cold limit and 
rapidly decays at larger p. The momentum dependence of the order parameter for dipolar gases results in a nonuniform 
energy gap for single-particle excitations and can, for example, manifest itself in processes with a large (of the order 
oi Pf) momentum transfer. 




In Eq.(p3|) the anisotropy of the order parameter in the momentum space is described by the function cf)o{'n) = 
V^sin [(7r/2) cos^]. The order parameter is an odd function of cos9 and is negative for 7r/2 < 9 < n. This does not 
cause any problems as all physical quantities are determined by |A|. The maximum value of |A(0,n)| is reached in 
the direction of the dipoles, i.e. for 6* = (0o = V^)- In the direction perpendicular to the dipoles {6 — n/2) the 
order parameter vanishes. 

If we consider only the p — p scattering channel, the angular dependence of the order parameter will be determined 
by the function a/3 cos instead of (/)o(n). The coefficient Cqoo is then equal to 9/4, and the result for the order 
parameter reads: 

App(0,n) = l^yi^Tr^y^ZlVacosf? = 2.5r,y^^V2cos0; « 1- (34) 

The angular dependence of |App(0, n)| is qualitatively similar to that of the true order parameter |A(0,n)|. The 
maximum value of |App(0, n)| is also reached in the direction of the dipoles and it is exactly equal to the maximum 
value of |A(0, n)|. Also, |App(0, n)| vanishes in the perpendicular direction. However, for intermediate values of 9 the 
quantity |App(0, n)| can be up to 40% smaller than |A(0,n)|. 

The anisotropy of the order parameter provides a major difference of the properties of the superfluid dipolar Fermi 
gas from those of the (two-component) fermionic gas with the s-wave pairing due to short-range intercomponent 
interaction. This anisotropy ensures the anisotropic momentum dependence of the gap in the spectrum of single- 
particle excitations, which appears below the transition temperature Tc- For example, excitations with momenta in 
the direction of the dipoles acquire the largest gap. In contrast to this, the eigenenergies of excitations with momenta 
perpendicular to the dipoles remain unchanged. The properties of collective excitations are also expected to have a 
nontrivial dependence on the direction of their momenta. Therefore, the response of the dipolar superfluid Fermi gas 
to small external perturbations will have a pronounced anisotropic character. 

Another distinguished feature of the superfluid dipolar Fermi gas is related to the temperature dependence of the 
specific heat. Well below the critical temperature the single-particle contribution to the specific heat is proportional 
to T^, rather than being exponentially small as in the case of the s-wave pairing. This follows from the fact that 
the energy e of single-particle excitations has a line of zeros on the Fermi surface: e{pF) = for the angles at which 
A{pp,n) = 0, i.e. for 9 = n/2 and an arbitrary azimuthal angle ip. As a consequence, the density of states in the 
vicinity of the Fermi energy is z^(e) ~ e for e ^ Aq. Therefore, at temperatures T ^ T^, the temperature 

dependent part of the energy of the system is proportional to T^, and the specific heat is hence proportional to T^. 
This contribution is much larger than the one of collective modes which is cx and is dominant in the case of the 
s-wave pairing. 

It should also be mentioned that the properties of the superfluid dipolar fermionic gas are different from the 
properties of the gas with the p-wave pairing originating from a short-range attractive interaction in the p -wave 
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channel. The reason is that in the latter case the order parameter is isotropic, similarly to the -B-phase of superfluid 
•^He. The order parameter of dipolar gases is also different from the order parameter of the A phase of ^He where the 
gap vanishes only at two points on the Fermi sphere, i.e. at the poles of the sphere 9 = and 9 = tt. 

The anisotropy of the order parameter of dipolar gases is similar to that in the polar phase of superfluid liquid ^He, 
not realized in experiments as it has higher energy than experimentally observed A and B phases (see, e.g. Ref. p2t). 
For the polar phase the order parameter is also equal to zero on the equator of the Fermi sphere {9 = 7r/2 and an 
arbitrary ip) . The situation where the order parameter is zero on one or several lines on the Fermi surface is encountered 
in heavy-fermion compounds (for a review of possible superconducting phases of heavy-fermion compounds belonging 
to different crystalhne groups see, e.g. Ref. In these cases the temperature dependence of the specific heat is 

also (X (see, e.g. Ref. ^4]. However, from a general point of view, one would expect a different physical behavior 
of dipolar gases, for example with regard to the frequency and angular dependence of the response. This is because of 
the different types of symmetry groups: continuous rotational group for dipolar gases and discrete crystalline group 
in the case of heavy-fermion materials (see Ref. |^ for more details) . 

VII. CONCLUDING REMARKS 

Our results show prospects for achieving the BCS transition in single-component trapped gases of dipolar particles, 
in particular for (electrically polarized) fermionic polar molecules. As has been shown in Refs. |l^ , p5| , p6| , for trapped 
gases the BCS transition temperature is close to that of the 3D uniform gas with density n equal to the maximum 
density in the trap. This requires Tc to be much larger than the maximum trap frequency, which is generally the case 
for achievable temperatures. Therefore, we will use Eq.(p9|) for estimating Tc in the trapped case. 

We first compare our equations ( pT] ) and ( p9[ ) with the well-known BCS formula (see e.g. [|l8|) for the two-component 
Fermi gas with short-range attractive intercomponent interaction. In the latter case the exponent is expressed in terms 
of the s-wave scattering length a and is equal to T:h/2pF \a\. We then see that our dipole-dipole scattering with odd 
orbital angular momenta is equivalent to having the s-wave scattering length 

2m(f , , 

= (35) 



Accordingly, Eq. (p9|) takes the form 



1.44e;^exp<i--^^;>. (36) 
2pF \ad\ 



Qualitatively, the result of Eq.(|35|) is more or less expected, since |ad| turns out to be of the order of the characteristic 
radius of the dipole-dipole interaction, ~ mcP/h?', introduced in Section III. 

For most polar molecules the electric dipole moment ranges from 0.1 to 1 Debye. For example, the dipole moment 
of fermionic ammonia molecules ^^NDa |^ is d = 1.5 D, and we have the effective scattering length is = —1450 
A. This is larger than the scattering length for the intercomponent interaction in the widely discussed case of two- 
species fermionic gas of ^Li. From Eq.(^6|) we find that the BCS transition temperature for the single-component ND3 
dipolar gas will be larger than 100 nK at densities 71 > 5 • 10"'^^ cm~'^. Another interesting example is a linear fermionic 
molecule HCN which has dipole moment d = 2.98 D, and the corresponding effective scattering length aa = —7400 A. 

Remarkably, in ultra-cold single-component fermionic gases one can hope to reach much higher densities than in 
Bose gases and think of achieving the BCS transition for a smaller effective scattering length ad- The reason is that 
inelastic decay processes will be strongly suppressed due to the Pauli exclusion principle. For two identical fermions 
with momentum p of the relative motion, the pair correlation function behaves as (pr/Ti)^ at interparticle distances 
r smaller than the de Broglie wavelength h/p. Generally, inelastic processes occur at short interparticle distances Rq 
which in the ultra-cold limit are much smaller than h/p. Therefore, two-body inelastic collisions will be suppressed 
as (pRo)^ compared to the bosonic case where the pair correlation function is of order unity at any r outside the 
region of interparticle interaction. As a result, in a non-degenerate gas of fermions the inelastic rate decreases with 
temperature and is suppressed as T/eq, where the energy Eq = h^/mRQ. In a quantum degenerate Fermi gas a 
characteristic momentum of particles is of the order of pi? and the suppression factor is ~ {ep/so)- The suppression 



^ Bosonic (^^NDs) and fermionic (^^NDs) ammonia molecules were recently trapped at temperatures of around 30 mK in the 
Rijnhuizen experiment [Q. 
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of two-body inelastic collisions in fermionic gases was first found for spin relaxation in atomic deuterium [|j . For the 
rate of 3-body recombination we expect even a stronger suppression, i.e. by a factor of (T/eo)^ for the non-degenerate 
gas, and by a factor of [ep/e^Y the regime of quantum degeneracy (T < sp)- The suppression factor for inelastic 
collisions in single-component quantum degenerate Fermi gases, (eF/eo), is of the order of {nR^Y/^. The distance 
i?o is commonly smaller than 50 A. Therefore, even at densities n ~ 10"'^^ cm~^ one expects the suppression of two- 
body inelastic rates at least by 2 orders of magnitude, and the suppression of 3-body recombination by 4 orders of 
magnitude. This allows us to think of achieving densities n ^ 10^^ cm~'^ or somewhat larger, which is by more than 
a factor of 10 higher than the densities currently reached with ultra-cold Bose gases. 

For n ^ 10^^ cm~'^, the BCS transition temperature Tc is in the nanokelvin or microkelvin regime for dipolar 
fermionic gases with an effective scattering length ad ranging from —20 to —30 A. Such values of ad one easily finds 
in molecular gases. For example, fermionic ^^N^^O molecules have dipole moment d — 0.16 D and Eq.(p5|) gives 
Qd « -24 A. 

Interestingly, the effective scattering length \ad\ of the order of several tens of angstroms can be obtained in gases 
of atoms with induced dipole moments |2^j2^. By using a high dc electric field 10^ V/cm) one can induce 
permanent atomic dipole moments close to O.ID. The same values of d one obtains for the time averaged dipole 
moment of an atom, induced by a stroboscopic laser coupling of the ground atomic state to a Rydberg state | |T^ . The 
corresponding scattering length ad can then be close to —30 A, and one can think of achieving the BCS transition in 
such single-component atomic dipolar gases at densities n ^ 10^^ cm"'^ and temperatures ~ 100 nK. 
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